Enhanced zero-bias Majorana peak in disordered multi-subband quantum wires 
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A recent experiment [Mourik et al, Science 336, 1003 (2012)] on InSb quantum wires provides 
possible evidence for the realization of a topological superconducting phase and the formation of 
Majorana bound states. Motivated by this experiment, we consider the signature of Majorana bound 
states in the differential tunneling conductance of multi-subband wires. We show that the weight 
of the Majorana-induced zero-bias peak is strongly enhanced by mixing of subbands, when disorder 
is added to the end of the quantum wire. We also consider how the topological phase transition is 
reflected in the gap structure of the current-voltage characteristic. 
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Introduction. — A recent experiment [1 reports the re- 
alization of proximity-induced topological superconduc- 
tivity (2j-|4] and the formation of Majorana bound states 
in InSb quantum wires. Following theoretical suggestions 
|B] , superconducting order is induced in an InSb quan- 
tum wire by proximity to a Nb lead attached alongside 
the wire. At the other end, the quantum wire is contacted 
to a normal lead via a gate-induced tunnel junction. Ev- 
idence for the formation of Majorana bound states is 
found through measurements of the differential conduc- 
tance, which exhibits a zero-bias peak when a magnetic 
field is applied in certain directions. Similar results were 
also obtained for normal-metal-superconductor struc- 
tures based on InAs quantum wires |7_ . 

At zero temperature and in single-subband quantum 
wires, the Majorana-induced zero-bias peak is predicted 
to have a height of 2e 2 /h [HI [9] • At finite temperature, the 
zero-bias peak broadens with its weight fixed, so that the 
peak height is no longer expected to reach 2e 2 /h. In this 
paper, we consider the current-voltage characteristic of 
multi-subband wires — a situation which is presumably 
relevant to the experiment of Ref . pQ — , with a particu- 
lar emphasis on the dependence of the zero-bias peak on 
subband mixing by disorder. We show that, remarkably, 
the weight of the Majorana-induced zero-bias peak is typ- 
ically enhanced as the tunnel junction becomes more dis- 
ordered. The basic idea is that disorder couples the topo- 
logical channel, which itself is only weakly transmitted 
through the barrier, with the other non-topological sub- 
bands which have higher transmission coefficients. This 
coupling broadens the conductance peak and hence, in 
the presence of a finite temperature, enhances the zero- 
bias conductance [TU] . The intentional inclusion of disor- 
der in or near the barrier, either during the fabrication 
process of the InSb nanowires used in the experiment, or 
after fabrication of the device, could thus lead to an ad- 
ditional, strong signature of the Majorana end state. A 
similar effect is expected if the tunnel barrier is replaced 
by a point contact [11 , provided the point contact is 
non-adiabatic. 

The interest in Majorana bound states in low- 



dimensional condensed matter systems [TJ |7J H2TTTT] is 
driven by their remarkable properties: They are their 
own antiparticle, have zero energy, and obey non-Abelian 
exchange statistics [18, 19 upon adiabatic permutation of 
their positions. The latter two properties make Majorana 
bound states potentially useful for topological quantum 
computation |20j . 

Model system. — We consider a geometry close to that 
of the experiment in Ref. [T] shown schematically in 
Fig. [iji. It consists of a two-dimensional multi-subband 
semiconducting wire with spin-orbit velocity a, chemical 
potential /i, and width W. At one end, the semiconduc- 
tor is coupled laterally to a superconducting lead. At the 
other end, it is contacted to a normal metal via a tun- 
nel barrier defined by the gate potential U. The system 
is placed in a magnetic field parallel to the wire direc- 
tion with Zeeman energy B. Taking the x direction to 
be along the wire, the system is then described by the 
Bogoliubov-de Gennes Hamiltonian [H H21 ET| 




Figure 1. (Color online) (a) Setup of multi-subband quan- 
tum nanowire (NW) with gate-induced tunnel barrier (G) and 
proximity coupled s-wave superconductor (S). As in Ref. [T] 
we consider the conductance between the normal lead (N) 
and the superconductor. Subband mixing is induced through 
disorder in the short segment of length L between the tunnel 
barrier and superconductor, (b) Normal-state dispersion in 
the absence of disorder for four subbands with B = 1 meV 
and ma 2 /2 = 50 /ieV. 
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where the Pauli matrices a and r operate on the spin 
and particle-hole degrees of freedom, respectively. The 
parameter a y is included for future reference and equals 
a for the case of Rashba spin-orbit coupling. The lateral 
contact to the superconductor covers the region x > 0, so 
that we set A(x) = A for x > and A(x) = 0, otherwise, 
where A is the proximity-induced gap for B = 0. The 
disorder potential Vdis( r ) is nonzero in the region — L < 
x < between the gate-defined tunnel barrier and the 
superconducting contact only. In this region, we choose 
a Gaussian random potential with (Vdi s (r)) = and 

(V dis (r)V dis (r')) = ^S(r-r'), (2) 

where vf = y/2fjt/rn and fad are the Fermi velocity and 
mean free path. It is important to note that the mech- 
anism discussed here is different from reflectionless tun- 
neling [22 induced by disorder on the normal side of NS 
tunnel junctions. 

We numerically calculate the normal and Andreev re- 
flection matrices r oc (e) and rh c (e) for the Hamiltonian 
0, using the technique described in Ref. [23]. The dif- 
ferential conductance G(V) is then evaluated according 

to pi] 

G(V) = ir [1 + r hc (eV)r hc (eVy - r ee {eV)r ee (eV)^], 
h 

where the trace is in spin and channel space. In a wire 
of width W lateral momenta are quantized as p VjTl = 
Hmr/W with n = 1, 2, . . .. In our numerical calculations 
we use an effective mass m — 0.015m e , m e being the 
bare electron mass, proximity induced gap A = 250 /ieV, 
spin-orbit energy ma 2 /2 = 50 /J.eV, and width W = 110 
nm. [55]. This choice corresponds to the parameters of 
the InSb quantum wires used in Ref. pQ. The chemical 
potential in the nanowire is chosen as /i = 32.1 meV, 
corresponding to N = 4 occupied channels (cf. Fig. [lja) . 
The subbands in the nanowire are therefore separated by 
several meV and for Zeeman energies less than 1.5 meV 
(corresponding to B < 1 T for InSb) only the highest 
channel (subband index n = N = 4) can be in the topo- 
logical phase. 

Clean multi- subband quantum wires. — We first con- 
sider a clean multi-subband wire with a y — 0. To a good 
approximation, a gate-induced tunnel barrier exposes the 
electrons to a potential which depends only on the coor- 
dinate x along the wire. Consequently, the tunnel barrier 
does not mix the transverse subbands (channels) of the 
quantum wire and the subbands can effectively be con- 
sidered as independent. Each subband is characterized 
by a Fermi velocity Vp,n = ^(2m/j -p^) 1 / 2 , an exci- 
tation gap A„, and a transmission coefficient T n of the 
gate-induced tunnel barrier. Since the highest occupied 
subband n = N determines whether the wire is in the 
topological phase, we refer to this subband as the "topo- 
logical subband." A nontrivial topological phase exists if 
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Figure 2. (Color online) (a) Zero-bias conductance peak at 
zero temperature in a quantum wire with B — 0.5 meV and 
N = 4 transverse subbands, one of which is in the topological 
phase with barrier transmission T4 = 0.01. The three non- 
topological subbands have transmissions 2OT4, IOT4, and 4T4. 
The red curves show the conductance for four different disor- 
der configurations with I — 10L. The black dashed line shows 
the peak shape for the clean wire, (b) Same as in (a), but for 
a temperature T = 60 mK, larger than the zero-temperature 
peak width. 

B 2 > B 2 N = A 2 + (fi~pl_ N /2m) 2 016]. At the topologi- 
cal phase transition, i.e., for B — Bjq , the topological gap 
Ajy vanishes, whereas the excitation gaps for the other 
subbands remain finite. 

The topological phase is characterized by a zero-bias 
conductance peak 

with width T = /SoAjvTjv if L <C £ = hvp,N I Ajv • The 
numerical constant /3o takes the value /3q ~ 0.375 for 
the range of parameters we investigated (Tjy -C 1, A 
between and 60 /ieV). This width may be very small, 
since the transmission coefficient of the topological 
subband is typically much smaller than the transmission 
coefficients of the other channels. (For the iVth subband 
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Figure 3. (Color online) Probability distribution of the zero- 
bias peak width F in the presence of disorder in the nanowire 
segment — L < x < between the superconducting part and 
the barrier for a multi-channel wire with the same choice of 
parameters as in Fig. [2] With increasing disorder, T increases 
on average (red and green curve) due to subband mixing. For 
L 3> I Anderson localization reduces the overall transparency 
of the junction, causing F to decrease again in the case of very 
strong disorder (blue curve). 



to be topological, it is important that its band bottom be 
close to the chemical potential.) At finite temperature, 
the conductance peak is thermally broadened, 



G(V,T) 



deG(e,0)^f(eV-e,T), (4) 
as 



where /(e, T) is the Fermi distribution function at tem- 
perature T. Thermal broadening preserves the weight of 
the zero-bias peak. For k-gT ^> T, the peak width is of 
order fc B T, whereas the height (2e 2 /h)(TrT/4:k B T) is in- 
versely proportional to temperature. Both regimes are 
illustrated in Fig. [2] 

Effect of disorder. — If the disorder is limited to the 
segment of the semiconductor wire that is not in contact 
with the superconductor, i.e., to — L < x < (cf. Fig. [TJ , 
it has no effect on the existence of the topological phase 
[26) . However, impurity scattering in the "normal" part 
of the wire has profound consequences for the weight of 
the zero bias peak associated with the existence of the 
topological phase. The underlying reason is the large 
disparity in the transparencies of the different subbands, 
with the topological subband having the smallest trans- 
parency Tjy. Mixing of subbands by impurity scattering 
allows for the coupling of the topological subband to the 
normal lead via the lower subbands with higher trans- 
parency. 

The effect is illustrated in Fig. [2] where we show the 
shape of the zero-bias peak for various disorder configu- 
rations, such that the distance L between gate-induced 
tunnel barrier and the superconducting contact equals 



Figure 4. (Color online) Ensemble average of the contribu- 
tion Ti from disorder-induced subband mixing to the width 
r of the zero-bias peak as a function of disorder strength in 
the segment — L < x < 0. The peak width is normalized by 
the normal state conductance Gb = (2e 2 /Ii)NTb to focus on 
the effects of subband mixing and to eliminate changes in the 
overall transparency by Anderson localization. Inset: Contri- 
bution F2 to the peak width from lateral spin-orbit coupling 
for a rectangular barrier (red crosses) and a Gaussian bar- 
rier (blue dashed line). In both figures the parameters of the 
barrier potential have been chosen such that only T± differs 
appreciably from zero. 



one tenth of the characteristic scattering length 



I = hdVF,lVF,N/v 



(5) 



A systematic dependence on disorder strength can be 
seen in Figs. [3] and [4] Figure [3] shows the probability 
distribution of the zero-temperature peak width for dif- 
ferent values of the ratio L/l. We conclude that already 
a moderate amount of disorder causes subband mixing 
and an increase in the peak width. At very strong dis- 
order, L ^> I, Anderson localization suppresses the over- 
all coupling to the normal lead, leading to a decrease of 
the weight of the zero-bias peak. This effect is not re- 
lated to subband mixing and can be removed by normal- 
izing the peak weight to the normal-state conductance 
G B = (2e 2 /h)NT B of the device; see Fig. g The av- 
erage peak width from disorder-induced subband mixing 
saturates for L/l ^ 1 to (Ti) = PxTg/S.^, where j3\ 
is a numerical factor of the order of 0.1, the exact value 
depending on the barrier transparencies for different sub- 
bands and spin mixing due to the magnetic field. 

Other causes of subband mixing. — The lateral spin- 
orbit term proportional to a y in Eq. ([IJ may be an ad- 
ditional source of subband mixing. For small a y its con- 
tribution to the width of the zero-bias peak is 

r 2 = ftCm^, C = W ^ TB » ( 6 ) 
" W z kpVp^M 

proportional to with a numerical prefactor (3% that 
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depends on the precise shape of the barrier. In the in- 
set of Fig. [4] we show its effect on the conductance of 
a clean wire for a long and low tunnel barrier, so that 
only the lowest subband n = 1 has an appreciable trans- 
mission. For a rectangular barrier, the subband mix- 
ing caused by lateral spin-orbit coupling is maximal, but 
still weak in comparison to the maximal subband mix- 
ing obtained from disorder, since for B = 0.5 meV we 
obtain C ~ 10~ 4 . For a smooth barrier, which is the 
experimentally relevant limit, the numerical prefactor /3 2 
becomes vanishingly small and lateral spin-orbit coupling 
does not give any appreciable subband mixing. Subbands 
may also be mixed by a gate-defined barrier that is not 
perpendicular to the direction of the wire. The mixing 
effect is maximal if the barrier is rectangular, and effec- 
tively absent for smooth barriers. 

Current-voltage characteristic and topological gap. — 
Unlike in single-channel models for spinless p-wave su- 
perconductors, multi-subband models are characterized 
by the coexistence of multiple superconducting gaps in 
different sections of the Fermi surface. Specifically, the 
proximity-induced gaps in the lower subbands are only 
weakly affected by the applied magnetic field. In con- 
trast, the highest occupied subband should have a gap 
closing when it enters into the topological superconduct- 
ing phase at the critical magnetic field. Thus, it is in- 
teresting to investigate to which degree the differential 
conductance contains signatures of the gap closing at the 
topological phase transition and how disorder near the 
barrier affects these signatures. 

In Fig. [5^,, we show the differential conductance ver- 
sus bias voltage for a clean multichannel quantum wire 
at T = 60 mK. At the critical Zeeman field of the top- 
most channel B c = 0.27 meV a peak appears at zero 
bias voltage. Since the topological channel is only weakly 
transmitted through the barrier, its contribution to the 
conductance is weaker than that of the other channels. 
In conjunction with density of states effects [27], this ex- 
plains the very weak signature of the topological gap clos- 
ing in transport in Fig. [5^, consistent with the absence of 
the topological gap in the experimental measurements of 
Ref. pQ. As for the zero-bias peak, the gap-closing feature 
in the differential conductance will also be significantly 
enhanced by disorder in the barrier region. This is shown 
in Fig. [5J3 where both the zero-bias peak and the peaks 
associated with the topological gap for B < B c are much 
more pronounced than in Fig. [5^. Indeed, the topolog- 
ical gap originates from the same subband as the zero- 
bias peak and its visibility is thus enhanced by the same 
mechanism. Given that the predictions of the multiband 
model (JT|) are consistent with the experimental data of 
Ref. [J, the deliberate introduction of subband mixing 
would be an instructive probe of Majorana bound states. 

We acknowledge useful discussions with J. Alicea, Y. 
Oreg, G. Refael, and M. Wimmer. We are grateful for fi- 
nancial support by the Deutsche Forschungsgcmcinschaft 




bias voltage (mV) 

Figure 5. (Color online) (a) Differential conductance vs bias 
voltage in a clean multichannel nanowire for increasing B from 
to 0.5 meV (750 mT in InSb) in steps of 0.02 meV with the 
realistic parameters [T] a y = 0, T — 60 mK (feT = 5 meV), 
L — 10 nm, and Tjv = 0.01. The B > traces are offset 
vertically for clarity. The formation of a Majorana fermion 
is reflected in the emergence of a zero-bias peak. The corre- 
sponding closing of the topological gap is hardly discernible 
due to the low transparency of the topological channel. For 
B = there are coherence peaks at the proximity induced 
gap A = 0.25 meV. For larger Zeeman fields the bulk gap 
of the lower channels is decreased consistently with expecta- 
tions, (b) Same as (a) but with weak disorder in the region 
— L < x < adjacent to the barrier. All traces are calcu- 
lated for the same disorder configuration with a scattering 
length I = 10L. The zero-bias peak and the signature of the 
topological gap closing are considerably enhanced. 
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